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Recent results of baryonic B decays from Belle are reported. This study is done by a 78 fb _1 
data sample, consisting of 85.0 ± 0.5 million BB pairs, collected by the Belle detector at the 
KEKB asymmetric energy e + e~ (3.5 on 8 GeV) collider. The results reported here include 
the first observation of the two-body decay B° — + pA c , the first hyperonic decay B° — > p\w~ , 
and first observations of B + -> ppn+, B° -» ppK%, and B+ -> ppK*+ . 



1 Introduction 

The Belle collaboration recently reported the observation of B + — > ppK^^, which is the first 
known example of B meson decay to charmless final states containing baryons. The three-body 
decay rate is larger than the rate for two-body decays (such as B — > pp^), We continue this 
study with a larger data set and search for other related baryonic decay modes. In the Standard 
Model, these decays proceed via b — > u(c) tree and b — > s(d) penguin diagrams. They may 
be used to search for direct CP violation and test our theoretical understanding of rare decay 
processes involving baryons^. 

BelleSl is a general purpose detector operating at the KEKB asymmetric e + e collider. Data 
sample used here consists of 85.0 ± 0.5 million BB pairs or ~78 fb _1 data set collected on the 
T(45) resonance. 



2 Analysis Procedure 

The biggest challenge of observing rare B decay processes is to fish out a few signal events from 
a huge sample of background events. For example, after the trigger and hadronic event pre- 
selection, there are about 170 million e + e~ — > qq continuum events and 85 million BB events 
left for the 78 fb _1 data set. It is a tough job to reject all of the background. The following is 
a brief description of the procedure of a typical analysis. 



2.1 Signal identification 

Since the center-of-mass energy is set to match the T(4S) resonance and T(45) decays into a 
BB pair, we can use the following two kinematic variables to identify the reconstructed B meson 
candidates: the beam-energy constraint mass, Mb c = J -Ef eam — y\ , and the energy difference, 
AE = Eb — Ei, eam , where E^ eam , ps and Eb are the beam energy, the momentum and energy 
of the reconstructed B meson in the rest frame of T(45), respectively. 

"Throughout this report, inclusion of charge conjugate mode is always implied unless otherwise stated. 



The resolution of Mb c is about 3 MeV/c 2 which is due to the spread of the beam energy. 
Typically, the resolution of AE is about 10 MeV for final states with charged particles only. 

2.2 Background suppression 

The generic B decay is mainly via the b — > c transition which normally has more final state 
particles than those of the rare decay modes reported here, thus the background from generic B 
decays is much less than that from the continuum process. Similar B decay processes close to 
the target mode should be checked carefully one by one because they might feed into the signal 
region. 

For the continuum events, they have quite different event topology (more back-to-back or 
jet-like) than that of BB events (more spherical) in the T(4f?) frame. We can select some shape 
variables to form a Fisher discriminant in order to reject the continuum background. 

The probability density functions (PDFs) of the Fisher discriminant and other uncorrelated 
kinematic variables (e.g. the angle between the B flight direction and the beam direction), 
are combined to form signal (background) likelihood £g( BG ), and a cut is then applied on the 
likelihood ratio C1Z = £s/(£s + £bg) m order to fish out the signal. 

2.3 Yield determination 

After the optimization of selection cuts, events in the candidate M^ c — AE region are used for 
yield determination. This can be done either by un-binned likelihood fit or binned fit. The 
signal PDFs are normally a Gaussian function for M\, c and a double Gaussian for AE with 
parameters determined by MC simulation. The background PDFs are Argus function for M^ c 
and a straight line for AE with parameters determined by sideband events or by continuum MC 
simulation. If there is no evidence of signal events, one can use the fit results to estimate the 
expected background, and compare this with the observed number of events in the signal region 
in order to set the upper limit on the yield at the 90% confidence level^. 

2.4 Systematic check 

If one is lucky to observe a new decay mode, the comparison between signal MC simulation and 
data is necessary. However, most of the systematic studies are limited by the small statistics 
of the signal events. Therefore, control samples with large statistics and relevant to the study 
are checked in order to determine the systematic errors. For example, the systematic error due 
to the efficiency error of the proton identification requirements is studied with the A — > pn~ 
sample; the kaon identification is studied with the D* + — > D°7r + , D° — > K~ir + sample; the 
tracking efficiency is studied with the r\ — ► 77 and n — ► tt + -k~it sample, etc.. 

3 Recent results 

3.1 Results of B° -> pA c 

Although the four- and three-body baryonic B decays of B — > pA c 7r + 7r~ and B + — > pA c ir + are 
experimentally well observed^, there has no two-body mode being found. An effort is made to 
search for B° — > pA c with A c — > pK + -K~ . Fig. ^ shows the results. The measured branching 
fraction is B(B° -» pA c ) = (2.19±g;|| ± 0.32 ± 0.57) x 10" 5 , where the first and the second 
errors are statistical and systematic,respectively. The last error is due to the uncertainty in the 
branching fraction B{A C — ► pK + Tr~y*. This is the first ever observation of a two-body baryonic 
B decays. One interesting feature is that this two-body decay rate is about factor of 10 smaller 
than that of the related three-body decay. In contrast, the two- and three-body mesonic B 
decays are comparable. 
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Figure 1: Candidate B° — > pA c events: (a) scatter plot of AE versus Mbc, (b) AE distribution for Mbc > 
5.27GeV/c 2 , and (c) M bc distribution for \AE\ < 0.03 GeV. Solid curves indicate the fit results. 



3.2 Results of B° -> pAn' 
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Figure 2: The fitted yield divided by the bin size for B° -> pAir~ as a function of M p ^. The shaded distribution 
is from a phase-space MC simulation with area normalized to signal yield. 

We observe a clear signal for B° — > pA7r~, which is the first hyperonic 5 decay being 
found. Since the decay is not uniform in phase space, we fit the AE signal yield in bins of 
M p A, and correct for the MC-determined detection efficiency for each bin. This reduces the 
model dependence of the branching fraction determination. The signal yield as a function of pA 
mass is shown in Fig. [21 This threshold peaking behavior is similar to that of B + — > ppK + ^. 
The measured branching fraction is B(B° -> pAvr") = (3.97to° 8 o(stat.) ± 0.56(syst.)) x 1CT 6 . 
Searches for B — > pAK~ and pS°7r~ yield no significant signals and we set 90% confidence-level 
upper limits of B(B° -» pAlf - ) < 8.2 x and -» pS°vr-) < 3.8 x 1(T 6 . 

3.3 Preliminary results of B — > pph^ 

Following our observation of i? + — > ppK + mode, new theoretical models tried to explain the 
experimental data and illustrate the possible correlations between similar decay modes. Some 
theoretical explanation^! about the threshold peaking behavior suggests possible glue-ball in- 
duced enhancement. Table ^ summarizes the search results of related three-body decays. There 
is no evidence for possible glue-ball states in 2.2 GeV/c 2 < M p p < 2.4 GeV/c 2 region. These 
results are preliminary. Note that we apply charm veto^ to these decay modes. However, the 
vetoed events with pp coming from J ftp decay can be used as a calibration tool. The determined 
branching fractions using these events and the PDG 7 listed values are also shown in Tabled for 



comparison. 



Table 1: Summary of pph^ results. Branching fraction products, B(B — > J/tph^) x B(J/tp — > pp), are indicated 
with the J/V' symbol in the table. The search for B + — > glue — ball.R" + , glue — ball — > pp yields an upper limit at 

90% C.L. which is listed in the last row. 





£(xl(T6) 


Sig. 


J/^(xl0 - 6) 


PDG J/^(xl0~6) 


ppK + 


5.66+^ ± 0.62 


15.3cr 


2.48 ±0.25 


2.12 


Pp7T + 

PpKs 


3.06lg;^ ± 0.37 
0.94ta45 ±0.12 


6.7o- 
5.1cr 


1.25 ±0.25 


0.92 


ppK* + 


10.31±|;|| ± 1.55 


6.0cr 


2.41 ± 1.12 


2.94 


ppK*° 


< 7.6 x 10~ 6 UL at 90% C.L. 




2.97 ±0.57 


2.78 


glue-ball 


< 0.9 x 10" 6 UL at 90% C.L. 









Acknowledgments 

The author wish to thank the KEKB accelerator group for the excellent operation of the KEKB 
accelerator and the Moriond QCD organization committee for making such a wonderful confer- 
ence. This work is supported by the National Science Council of the Republic of China under 
the grant NSC-91-2112-M-002-028. 

References 

1. K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 88, 181803 (2002). 

2. K. Abe et al. (Belle Collaboration), Phys. Rev. D 65, 091103 (2002). 

3. For some of the recent developments, see: H.Y. Cheng and K.C. Yang, Phys. Rev. D 65, 
054028 (2002); H.Y. Cheng and K.C. Yang, Phys. Rev. D 66, 014020 (2002); C.K. Chua, 
W.S. Hou and S.Y. Tsai, Phys. Rev. D 66, 054004 (2002). 

4. A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A 479, 117 (2002). 

5. J. Conrad et al, Phys. Rev. D 67, 012002 (2003). 

6. SA. Dytman et al. (CLEO Collaboration), Phys. Rev. D 66, 091101 (2002); N. Gabyshev 
et al. (Belle Collaboration), Phys. Rev. D 66, 0911102 (2002). 

7. K. Hagiwara et al. (Particle Data Group), Phys. Rev. D 66, 010001 (2002). 

8. C.K. Chua, W.S. Hou and S.Y. Tsai, Phys. Lett. B544, 139 (2002). 



